Auxiliary Field Diffusion Monte Carlo calculation of nuclei with A < 40 with tensor 

interactions. 
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We calculate the ground-state energy of 4 He, 8 He, le O, and 40 Ca using the auxiliary field diffusion 
Monte Carlo method in the fixed phase approximation and the Argonne v' 6 interaction which includes 
a tensor force. Comparison of our light nuclei results to those of Green's function Monte Carlo 
calculations shows the accuracy of our method for both open and closed shell nuclei. We also apply 
it to 16 and 40 Ca to show that quantum Monte Carlo methods are now applicable to larger nuclei. 



The nuclear many-body problem has not found a sat- 
isfactory solution yet. As opposed to what happens for 
other systems, one has to face both the absence of a 
definitive scheme for the description of nucleon-nucleon 
forces, and with the extreme complexity induced by the 
fact that forces are state dependent. Recently there have 
been a few attempts to reduce the problem to a more 
fundamental level by exploiting in full the scheme of Ef- 
fective Field Theory (EFT). However, this approach, at 
present, is applicable only to very small nuclei In- 
tegration of subnuclear degrees of freedom in EFT now 
provide interactions which are almost as accurate in de- 
scribing scattering data as the most popular realistic 
interactions [2j- However, the assessment of the qual- 
ity of different potentials in many-body sistems remains 
problematic. Quantum Monte Carlo is the only available 
tool that provides estimates of physical observables with 
an accuracy comparable to that obtained by means of 
few-body techniques [1[ , but in a wider range of nuclear 
masses. QMC calculations might therefore help to dis- 
cuss and gauge the validity of the proposed interactions 
without the bias implied by the use of approximate meth- 
ods. However, the operatorial structure of such poten- 
tials, both phenomenological and EFT, prevented to push 
calculations beyond A = 12 QMC methods deal with 
the exponential increase in the computational time with 
particle number by sampling the degrees of freedom to 
evaluate sums and integrals. The variational and Green's 
function Monte Carlo methods [3] are some of the most 
successful methods for calculating the properties of light 
nuclei. There the spatial degrees of freedom are sampled, 
but the spin isospin degrees of freedom of the nucleons 
are explicitly summed and not sampled. The exponential 
growth of the spatial degrees of freedom is controlled, but 
since there are four spin-isospin states per nucleon, the 
computations grow exponentially - roughly as four raised 
to the number of nucleons. To make the method com- 
putationally efficient, the spin-isospin degrees of freedom 
must also be sampled. That is what the Auxiliary Field 



Diffusion Monte Carlo (AFDMC) [fa] does in the most ef- 
ficient way known today. 

Here we will demonstrate that AFDMC can be used to 
solve for the energy of nuclei with nucleons interacting 
via a somewhat simplified two-body interaction which, 
however, does contain the tensor interaction that comes 
from the one-pion exchange potential, and already in- 
cludes all the terms which make standard approaches 
unpractical. This interaction term is indeed the most 
demanding one, together with the short range repulsion, 
in solving the nuclear many-body Schrodinger equation. 
Its inclusion in the Hamiltonian provides one of the most 
severe tests of the efficiency of a quantum Monte Carlo 
algorithm. Adding the neglected spin-orbit terms and 
the three-body potential is not expected to change the 
main conclusion of this paper, namely that AFDMC is 
applicable with the same accuracy to both nuclei and 
nuclear matter. 

The AFDMC method which includes a path constraint 
to control the fermion sign problem, has given good re- 
sults for pure neutron matter @, 0], for neutron-dropsjsj 
and for the valence neutrons of neutron-rich nuclei[9(. 
For nuclei, the strong tensor force in the isosinglet chan- 
nel makes sampling the spin-isospin states more difficult, 
leading to unsatisfactory results when np and pp inter- 
actions are active. In this letter we demonstrate that 
AFDMC in the fixed phase approximation overcomes this 
problem, and, as a consequence, can be applied to cal- 
culate binding energy of large nuclei and nuclear matter 
with realistic interactions. Most of our results are with 
the Argonne v' 6 interaction a simplified version of the 
Argonne vi& potential ll| truncated to the first six op- 
erators and modified to describe the binding energy of 
deuteron. We calculate the energy of the alpha particle 
and the open shell nucleus 8 He and test the accuracy of 
our results by comparing them to those from GFMC[lO(. 
We apply the same algorithm to the study of binding en- 
ergy of 16 O first with the Argonne V14 interaction trun- 
cated to include only six operators to compare our results 
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with those of Cluster Variational Monte Carlo [12| and 
Fermi Hypernetted Chain in the Single Operator Chain 
approximation (FHNC/SOC). Finally, we calculate, the 
ground-state energy of the closed-shell nuclei 16 O and 
40 Ca using the Argonne v' e interaction. 

Monte Carlo methods are most efficient when applied 
to sums and integrals with positive kernels. These can 
be interpreted as probabilities and probability densities. 
We have chosen to use the Argonne series of potentials [Ioj| 
because they are substantially spatially local with only a 
few derivative terms. The short-time Green's functions 
that we need are then readily calculated and sampled. 

Our Hamiltonian is 
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where i and j label the two nucleons, Ti 3 is the dis- 
tance separating the two nucleons, and the in- 
clude spin and isospin operators, where M is the num- 
ber of operators (i.e. 18 in v±s models). The mass 
to -1 = (nip 1 +m~ 1 )/2 where m p and m n are the proton 
and neutron masses. We use the Argonne v' R model 
where the two-body potential is reprojected from the Ar- 
gonne vis to the M = 6 level. The six 0^ n '(i,j) terms 
are the 1 , n-fj, a t ■ a 3 , (<?; • a 3 ) (fj -Tj), S i3 , and Syr • t 3 , 
where Si 3 is the tensor operator 3(7^ • Ti 3 o 3 ■ fij — di ■ S 3 . 
The fi and di are the Pauli matrices for the isospin and 
spin of particle i. The inclusion of neutron- proton mass 
difference, electromagnetic interactions, spin-orbit inter- 
actions, and three-body potentials can be done with an 
increase in complexity. However the bulk of the bind- 
ing energy comes from the vq terms which include the 
one-pion exchange parts of the potential. 

Traditionally, ground-state quantum Monte Carlo cal- 
culations begin by using a variational calculation to op- 
timize a trial wave function. This trial wave function 
is then used to guide the sampling of the random walk 
in diffusion or Green's function Monte Carlo. A typi- 
cal form for a good trial function has a model function 
| $) given by a small linear combination of antisymmetric 
products of orbitals multiplied by a symmetrized product 
of two-body operator correlations. Evaluating this sym- 
metrized product trial function at the spatial positions 
R and spin-isospin values S gives the expression 
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Unfortunately, the evaluation of this wave function re- 
quires exponentially increasing computational time with 
the number of particles. Since the evaluation for all spin- 
isospin values has the same computational complexity, 
light nuclei variational and Green's function Monte Carlo 
calculations sum the spin-isospin degrees of freedom. 



Since for large numbers of particles we cannot evaluate 
these trial functions, we use much simpler wave functions 
which contain only the central Jastrow correlation. The 
evaluation of our simpler wave function require order A 3 
operations to evaluate the Slater determinants and A 2 
operations for the central Jastrow. Since many important 
correlations arc neglected in these simplified functions, 
we use the Hamiltonian itself to define the spin sampling. 

Specifically our trial function is 
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where A is an antisymmetrization operator, <p are single 
particle space and spin-isospin orbitals, built from com- 
binations of radial functions, spherical harmonics and 
spinors. i? cm = A^ 1 2^iLi ^ is the center of mass of 
the nucleus. The Jastrow function is the central part of 
the FHNC/SOC correlation operator Fi 3 which min- 
imized the energy of nuclear matter at /?o=0.16 fm . 

Radial orbitals are calculated in the self-consistent po- 
tential generated by the Hartree-Fock algorithm with 
Skyrme's effective interactions of Ref. [14( that has been 
used to study light nuclei. Given a set of positions and 
spinors, the antisymmetrization produces a determinant 
of single particle orbitals. For open-shell nuclei, a sum 
of several determinants is used to build a trial wavefunc- 
tion of a total angular momentum J that describes the 
nucleus. The determinants are multiplied by the central 
Jastrow factor to give the value of our trial function. 

The AFDMC method works much like diffusion Monte 
Carlo Hi, i, 0. The wave function is defined by a set 
of what we call walkers. Each walker is a set of the 3A 
coordinates of the particles plus A normalized four com- 
ponent spinors representing the spin-isospin state. The 
imaginary time propagator for the kinetic energy and the 
spin-independent part of the potential is identical to that 
used in standard diffusion Monte Carlo. The new posi- 
tions are sampled from a drifted Gaussian with a weight 
factor for branching given by the local energy of these 
components. Since these parts do not change the spin 
state, the spinors will be unchanged by these parts of 
propagator. 

To sample the spinors we first use a Hubbard 
Stratonovich transformation to write the propagator as 
an integral over auxiliary fields of a separated product 
of single particle spin-isospin operators. We then sample 
the auxiliary field value, and the resulting sample inde- 
pendently changes each spinor for each particle in the 
sample, giving a new sampled walker. 

Specifically we write the vq interaction as 
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Where v\ (r) is the central interaction, and V s d is the spin- 
isospin dependent part. The A matrices depend only on 
the positions of the particles. They are zero when i = j 
and they are real and symmetric so that they have real 
eigenvalues A„ ' , a! <TT ' ) , A„ and corresponding real nor- 
malized eigenvectors ipn(i,a), ip„ (i,ct),ipn (*)• The 
spin-dependent potential can be written as a sum of 
squares of single-particle operators as 
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where the 15^4 operators 
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are labeled in a convenient order of the eigenvectors and 
the corresponding eigenvalues. 

We apply the Hubbard-Stratonovich transformation to 
write 
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where x is an auxiliary field. Each of the 15A terms in 
Eq. [^requires an auxiliary field. We write the short time 
approximation of the spin-dependent propagator as 



-V sd At 



dX exp 



where dX = Yll, 5 ^ -^£s= and we drop commutator terms 
which are higher order than At on the right hand side. 

Once the Hubbard-Stratonovich variables have been 
sampled, the resulting propagator acting on a walker (i.e. 
positions and spinors) gives a single new walker. 

Since the trial function evaluated at the walker spin- 
isospin and position can be complex, we use a fixed-phase 
approximation [lij . We importance sample the auxiliary 
field variables x m by writing 
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where (O n ) = ($ T \O n \R, S)/($ T \R, S) is the mixed ex- 
pectation value. The first two terms are then combined 



TABLE I: The Ground-State energies of the alpha particle 
and of 8 He calculated with different methods using the Ar- 
gonne v' 6 interaction. The GFMC results are taken from Ref. 
[Tod after the subtraction of the Coulomb term of 0.7MeV[Tj§]. 
The EIHH result doesn't contain the Coulomb interaction. 
All the value are expressed in MeV. 



method 


£( 4 He) 


£( B He) 


AFDMC 
GFMC 
EIHH 


-27.13(10) 
-26.93(1) 
-26.85(2) 


-23.6(5) 
-23.6(1) 



to form a shifted contour Gaussian as in Ref. |l6|. Ap- 
plying the fixed phase approximation, instead of the pre- 
viously used constrained path, the walker weight can be 
reexpressed in terms of the local energy 



E L (R, S) = Re 



(V T \H\RS) 
(Vt\RS) 
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This change to the original algorithm was necessary to 
overcome the unphysical discrepancies observed in earlier 
AFDMC work on nuclei with tensor forces [l7j |. 

Our algorithm becomes: i) sample \R,S) initial walk- 
ers from | (^t\R, S) 1 2 using Metropolis Monte Carlo; ii) 
propagate in the usual DMC way with a drifted Gaussian 
for a time step; iii) diagonalize, for each walker, the po- 
tential matrices and A( crr ); iv) sample the cor- 
responding shifted contour auxiliary field variables and 
update the spinors. The new walker has a weight given 
byexp(-E L {B?,S')At). 

Our trial function contains no tensor correlations and 
the variational estimate is not even bound. The diffu- 
sion process enforced by the AFDMC method is capable 
of crossing the transition from an unbound to a bound 
system, leading to energy estimates which compare very 
well with the available GFMC results. Table U reports 
results for the alpha particle and the open shell nucleus 
8 He. For the alpha particle AFDMC estimates are com- 
pared with GFMC and the Effective Interaction Hyper- 
spherical Harmonic (EIHH) methods^. The AFDMC 
agreement with GFMC and EIHH for 4 He is within about 
1% of the total energy. The agreement between AFDMC 
and GFMC for 8 He is even better. 

We have compared our results for 16 O with other meth- 
ods. The variational FHNC/SOC[2l[, and Cluster varia- 
tional calculations [12| used the Argonne U14 interaction. 
Our result for the energy, keeping just the first six op- 
erators, is -90.8(1) MeV. The variational results keeping 
just those same six operators are -83.2 MeV from Vari- 
ational Monte Carlo and -84.0 MeV from FHNC/SOC. 
The AFDMC method seems to lower the energy by about 
10% with respect to the two different variational results; 
however, the variational wavefunctions were optimized 
with the full vu interaction instead of our truncated in- 
teraction. 
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TABLE II: Computed Ground-State Energy in MeV of 4 He, 
8 He, 16 and 40 Ca for the Argonne v'q interaction. Experi- 
mental energies are also reported [2^]. We also calculated the 
energy of 28 nucleons in a periodic box to extrapolated the 
nuclear matter energy at equilibrium density as described in 
Ref. [53] All the value are expressed in MeV. 



nucleus 


E 


E/A 


E eX p 


Eexp / A 


4 He 


-27.20(5) 


-6.8 


-28.296 


-7.074 


8 He 


-23.6(5) 


-2.95 


-31.408 


-3.926 


16 


-100.7(4) 


-6.29 


-127.619 


-7.98 


40 Ca 


-272(2) 


-6.8 


-342.051 


-8.55 


nuclear matter 




-12.8(1) 







We then performed calculations for the 16 O and for 
40 Ca with the Argonne v' 6 NN interaction. We chose this 
potential because it is a reprojected version of the more 
sophisticated Argonne vi$ . Its main deficiency is the lack 
of a spin-orbit interaction and the three-body potential. 

Results arc reported in table [Til where it is reported 
also the energy of nuclear matter at the equilibrium den- 
sity p =0.16fm- 3 calculated with AFDMC[H. 

As expected the v' 6 interactions is not sufficient to build 
the total binding energy of 16 O and of 40 Ca. This NN 
interaction gives about 96% of total binding energy for 
alpha particle, 75% for 8 He, 79% for 16 and 79% for 
40 Ca. Our 16 O is unstable to break up into 4 alpha par- 
ticles, and our 40 Ca has the same energy of 10 alpha 
particles. This is behavior is consistent with the simple 
pair counting argument of Ref. [24j . The surface energy 
coefficient in the Weizsacker formula, resulting from the 
comparison of the binding energies per nucleon of sym- 
metrical nuclear matter and 40 Ca is 20.5 MeV, not too 
far from the experimental value of 18.6 MeV. 

We have extended the AFDMC method, used to per- 
form calculations for neutron systems, to study finite nu- 
cleonic systems and have obtained agreement with results 
from other methods for light nuclei. The method eas- 
ily accommodates open-shell nuclei. We have calculated 
systems with up to A=40 nucleons here. The computa- 
tional time per imaginary time step scales as A 3 . The 
total computational time depends on the desired quan- 
tity, the distribution of excited states, and the quality 
and complexity of the trial function just as in all other 
quantum Monte Carlo methods. We believe that the re- 
sults we have presented show that AFDMC in the fixed 
phase approximation has become a very powerful tool 
to solve large nuclear systems with realistic interactions. 
This opens up the possibility of calculating at an enriched 
accuracy heavy nuclei and asymmetric nuclear matter. 

To be able to predict accurately the structure of nuclei, 
we must use a more realistic Hamiltonian. The main 
two features missing from this work are the three-body 
potential and spin orbit terms. 



The Urbana-IX potential has been used by us in 
previous neutron studies. This potential contains a 
spin-independent short range repulsion, and operator 
terms based on the Fujita-Miyazawa model. The spin- 
independent part and the so-called anticommutator 
terms lead to terms with two- or fewer spin-isospin opera- 
tors and can be included immediately in AFDMC calcula- 
tions. The commutator terms have as well as the isospin 
exchange spin-orbit require additional auxiliary fields to 
rewrite them in terms of single-particle spin-isospin oper- 
ators. The other terms in a more realistic interaction are 
included perturbatively in GFMC calculations. Initially, 
we will also try to include them perturbatively. This will 
likely be more difficult in AFDMC than in GFMC since 
our trial wave functions are much simpler and much less 
accurate so that the perturbative estimates will likely 
also be less accurate. The inclusion of these terms is in 
progress and will be the subject of future works. 

We thank G. Orlandini, W. Leidemann, and S.C. 
Pieper for providing us the EIHH and GFMC results with 
v' e for comparison. This work was supported in part by 
NSF grant PHY-0456609. Calculations were performed 
on the HPC facility "BEN" at ECT* in Trento under a 
grant for Supercomputing Projects. 
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